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ABSTRACT 


The Sloan Digital Sky Survey-Ill (SDSS-III) Apache Point Observatory 
Galactic Evolution Experiment (APOGEE) has obtained high resolution (R ~ 
22,500), high signal-to-noise ratio (> 100) spectra in the H—band (~1.5—1.7 
/xm) for about 146,000 stars in the Milky Way galaxy. We have computed 
spectral libraries with effective temperature (Tefr) ranging from 3500 to 8000 
K for the automated chemical analysis of the survey data. The libraries, used 
to derive stellar parameters and abundances from the APOGEE spectra in the 
SDSS-III data release 12 (DR12), are based on ATLAS9 model atmospheres and 
the ASSeT spectral synthesis code. We present a second set of libraries based 
on MARGS model atmospheres and the spectral synthesis code Turbospectrum. 
The ATLASO/ASSeT (Tefr = 3500—8000 K) and MARGS/Turbospectrum (Tefr 
= 3500—5500 K) grids cover a wide range of metallicity (—2.5 < [M/H] < +0.5 
dex), surface gravity (0 < log (? < 5 dex), microturbulence (0.5 < ^ 8 km s“^), 

carbon (—1 < [G/M] < +1 dex), nitrogen (—1 < [N/M] < +1 dex), and a-element 
(—1 < [a/M] < +1 dex) variations, having thus seven dimensions. We compare 
the ATLAS9/ASSeT and MARGS/Turbospectrum libraries and apply both of 
them to the analysis of the observed H—band spectra of the Sun and the K2 
giant Arcturus, as well as to a selected sample of well-known giant stars observed 
at very high-resolution. The new APOGEE libraries are publicly available and 
can be employed for chemical studies in the H—band using other high-resolution 
spectrographs. 


Subject headings: astrochemistry - radiative transfer - stars: atmospheres - surveys 



1. Introduction 


Previous near-infrared (JHK bands) spectroscopic observations of individual or small 
selected samples of giant stars have been limited in scope and mostly biased towards the 
brightest sources. This has prevented the study of the chemical abundance patterns in 
unbiased samples of stars towards the inner (and dusty) parts (e.g., Galactic bulge and 
center) of our Galaxy, signihcantly hampering progress towards a full understanding of the 
formation and chemical (and dynamical) evolution of the Milky Way. This unfortunate 
situation has dramatically changed in the new era of massive high-resolution spectroscopic 
surveys. In particular, the Apache Point Observatory Galactic Evolution Experiment 
(APOGEE) has focused on collecting high-resolution and high-quality H—band spectra for 
a large (> 10^ stars) sample of giant stars, with access to the inner - and more extinguished 
- regions of our Galaxy. 

APOGEE is one of the four spectroscopic surveys of the Sloan Digital Sky Survey 
III (SDSS-III; e.g., Eisenstein et ah 2011). It is a high-resolution (R = A/AA = 22,500) 
H—band (1.514-1.696 /xm) spectroscopic survey spanning all stellar populations in our 
Galaxy (see e.g., Allende Prieto et al. 2008; Eisenstein et al. 2011; Majewski et al. 2015). 
During the period from 2011 to 2014, the APOGEE survey collected about 500,000 spectra 
of ~146,000 stars, predominantly post-main sequence stars (red giants, subgiants, and 
red clump stars), using the Sloan Foundation 2.5m telescope (Gunn et al. 2006) and 
an innovative multi-object IR spectrograph (Wilson et al. 2010). High-resolution stellar 
spectra of red giants in the H—band show a rich diversity of absorption lines from a wide 
variety of atoms and molecules, with OH, GN, and GO the most important molecular 
contributors. To ascertain the stellar atmospheric parameters and measure chemical 
abundances from the observed spectra, the APOGEE Atmospheric Stellar Parameters and 
Ghemical Abundances Pipeline (ASPGAP; Garcia Perez et al. 2015) relies on an algorithm 



6 


that identifies the best-fitting synthetic spectrnm for each observed spectrum. The fitting 
code uses interpolation over pre-computed multi-dimensional grids of synthetic spectra (i.e., 
model stellar spectral libraries) to find the best model (with the minimum values) for 
each observed spectrum. Synthetic spectra are calculated using classical model atmospheres 
(see e.g., Meszaros et ah 2012) and extensive atomic and molecular line lists (Shetrone et 
al. 2015). 


The SDSS-III APOGEE public data release 10 (DRIO; Ahn et al. 2014) was based on 
synthetic spectral libraries computed using Castelli & Kurucz (2003) model atmospheres 
(see, e.g., Meszaros et al. 2013). The Castelli & Kurucz (2003) model atmospheres 
incorporate line opacity by means of opacity distribution functions (ODE) and are based 
on solar (or scaled solar) chemical abundances from Grevesse & Sauval (1998). For the 
APOGEE main survey targets, the fundamental stellar parameters (effective temperature 


(Teff), surface gravity (log g), and metallicity ([M/H]))I^ and the relative abundances of 
a-elements ([a/M]; in this case O, Mg, Si, S, Ca, and Ti), carbon ([C/M]), and nitrogen 
([N/M]) were released in DRIO. The final SDSS-III APOGEE public data release, DR12 
(Alam et al. 2015) is based on self-consistent spectral libraries, where the same chemical 
abundances are used both in the computation of the model atmospheres and in the spectral 
synthesis. Moreover, newer solar reference abundances from Asplund et al. (2005) are now 
adopted In addition to the main stellar parameters the individual element abundances 
for up to 15 elements (typically with a precision of 0.1 dex or better) are also released in 


^A linear relationship between microturbulence and surface gravity is adopted (see Garcia 
Perez et al. 2015); e.g., ^ = 2.478 — 0.325 x log g is used in DR12 (Alam et al. 2015). 

^The helium reference abundance adopted for the computation of the synthetic spectral 
libraries is: 12 -|- logig(NHe/NH) = 10.93, where Nne and Nh are the number density of 
helium and hydrogen nuclei, respectively. 
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DR12. 


In this paper we present for the first time the H—band stellar spectral libraries for the 
SDSS-III APOGEE survej^, which can be used as well for chemical studies using other 
high-resolution spectrographs working in the H—band. The DR12 spectral libraries are 
based on ATLAS9 model atmospheres (see Meszaros et ah 2012) and calculated with the 
ASSeT (Koesterke et ah 2008; Koesterke 2009) spectral synthesis code. In addition to 
the official family of ATLAS9/ASSeT DR12 spectral libraries, we have computed similar 
spectral libraries based on MARCS model atmospheres, with the Turbospectrum synthesis 
code (Alvarez & Plez 1998; Plez 2012). We provide a comparison between these different 
model atmospheres and spectral synthesis codes in order to check the validity of the adopted 
DR12 synthetic spectral libraries. In Section [2l we describe the parameter range and the 
calculation method of the ATLAS9/ASSeT and MARCS/Turbospectrum stellar spectral 
libraries, while Section 3 discusses systematic differences between the two grids of synthetic 
spectra. Both grids are applied to the observed H—band spectra of the Sun and the K2 
giant Arcturus in Section 4. In Section 5, we use both model stellar spectral libraries to 
derive the chemical patterns in a selected sample of well-known giant stars observed with 
the Fourier Transform Spectrograph (FTS) on the Kitt Peak National Observatory 4m 
Mayall reflector. Our main conclusions and future work are given in Section 6. 


2. Synthetic Spectral Libraries for APOGEE 

As mentioned above, the SDSS-III APOGEE DR12 results (Alam et ah 2015) are based 
on ATLAS9/ASSeT synthetic spectral libraries. The APOGEE synthetic spectral libraries 

^The stellar spectral libraries are available online; 

data.sdss3.org/sas/dr12/apogee/spectro/redux/speclib/asset/kurucz_filled/solarisotopes/ 




are continuously improved and will be updated in the future (see Section 6 ). In this section, 
we describe two H—band stellar spectral libraries developed for the SDSS-III APOGEE 
survey, based on ATLAS9 and MARCS model atmospheres and computed with the ASSeT 
and Turbospectrum spectral synthesis codes, respectively. All MARCS/Turbospectrum 
computations were performed on the Condor cluster at the Institute de Astrofisica de 
Canarias (lAC). Condor (or HTCondor) is a High Throughput Computing (HTC) system 
developed by the University of Wisconsin-Madison (UW-Madison)G. In the case of the 
lAC, Condor consists of a computer cluster with 808 CPUs. On the other hand, the 
ATLAS9/ASSeT calculations were performed on the clusters Stampede and Maverick 
operated by the Texas Advanced Computing Center (TACC). 


It is to be noted here that the ATLAS9 and MARCS model atmospheres do not include 
the line opacity for the polyatomic carbon molecules C 2 H 2 and C 3 (Meszaros et al. 2012). 
These molecules are known to dominate the infrared spectra of cool (Tgfj < 4000 K) carbon 
stars with C/0 > 1.0, strongly affecting their thermal atmospheric structure. Thus, at 
present the ATLAS9/ASSeT and MARCS/Turbospectrum synthetic spectra with Teg < 
3500—4000 K and C/0 > 1.0 are not reliable (see Section 6 ). 


2.1. ATLAS9/ASSeT spectral library 

The APOGEE ATLAS9/ASSeT spectral library makes use of the new ATLAS9 grid 
of model atmospheres presented in Meszaros et al. ( 2012 ). ATLAS9 model atmospheres 
(Kurucz 1993) are one-dimensional plane-parallel models computed in LTE and using 
opacity distribution functions (ODFs) to handle line opacity (see e.g., Kurucz 2005). The 
mixing-length scheme for convective energy transport is adopted, and the Kurucz atomic 

“^http: //research, cs.wise. edu/htcondor/ 
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and molecular line list^ are used (see Meszaros et al. 2012 for more details). 

The new ATLAS9 model atmosphere^ are based on the recent solar composition by 
Asplund et al. (2005). This has the advantage of matching the solar composition adopted 
also in the construction of the MARCS models described in the next Section. Meszaros 
et al. (2012) carried out new ODFs and Rosseland mean opacity calculations for several 
microturbulent velocities = 0, 1, 2, 4, and 8 km s“^), but all the APOGEE ATLAS9 
model atmospheres are computed with = 2 km s“^. 


For the purpose of building the libraries, we divided the ATLAS9 model atmospheres 
into two grids (GK and F) depending on the effective temperature (Tes), with the GK- 
and F-classes covering the 3500—6500 K and 5500—8000 K Tgfj ranges, respectively. 

It is important to note that missing (non-converged) ATLAS9 model atmospheres are 
substituted by the nearest model in chemical space at the same Tgff and log g. The 
non-converged ATLAS9 model atmospheres were a total of 6217 structures, most of them 
corresponding to cool high surface gravity models (Tgs < 4000 K and log > 4.0) and cool 
C-rich models (Tgfj < 4000 K and C/0 > 0.75 ). 

Most stars targeted by SDSS-III APOGEE are red giant and dwarf candidates, which 
makes the ATLAS9 GK-class model grid the most important for APOGEE. The ATLAS9 
GK-class models cover eleven and eleven log g values (from 3500 to 6000 K in steps 
of 250 K and from 0 to 5 dex in steps of 0.5 dex, respectively), seven metallicities [M/h|3 


®http://kurucz.harvard.edu/linelists.html 
®http://www.iac.es/proyecto/ATLAS-APOGEE/ 

^The overall metallicity [M/H] accounts for all elements with atom number Z > 2 and 
[M/H] = log]^Q(AM/A"H)*—logio(A^M/A"H) 0 , where Nm and Nh are the number density of any Z 
> 2 element and hydrogen nuclei, respectively. For the construction of the synthetic spectral 
libraries, all metals other than C, N, and a—elements are scaled with the corresponding 
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(from —2.5 to +0.5 dex in steps of 0.5 dex), and 81 values for the abundance of C and the 
a-elements ([C/M] and [a/M] from —1.0 to +1.0 dex in steps of 0.25 dex). The a-elements 
in the ATLAS9 models are as in the introduction. Thus, 68,607 GK-class ATLAS9 models 
are present in the GK ATLAS9 APOGEE grid used to produce the ATLAS9/ASSeT 
spectral library. 

The ATLAS9 F-class grid is used for the analysis of a much smaller number of warmer 
stars such as the telluric line standards observed by the SDSS-III APOGEE survey. The 
F-class ATLAS9-APOGEE grid contains also 68,607 model atmospheres that cover 
values from 5500 K to 8000 K (in steps of 250 K) and the same [M/H], [G/M], and [a/M] 
ranges as the GK-class one. 

The code ASSeT was used to compute synthetic spectra for the Kurucz ATLAS9 
models. ASSeT was originally written for calculating spectra for 3D hydrodynamical 
simulations, and later a ID branch was developed for dealing with large numbers of 
classical hydrostatic model atmospheres. The code has an option to pre-calculate the 
opacity as a function of temperature and a second thermodynamical quantity (density 
or electron density) on a grid covering the range of the models of interest, interpolating 
on-the-fly for solving the radiative transfer equation. Opacity interpolations were avoided 
in the calculations described above in order to maximize accuracy. ASSeT adopts the 
same software package used by Synspec/Tlusty (Hubeny & Lanz 1995; Hubeny 2006) for 
computing opacities, and the input data described by Allende Prieto (2008). Gontinuum 
opacities are mainly from the Opacity Project (Gunto et al. 1993) and the Iron Project 
(Nahar 1995; Bautista 1997), and line opacities are shared in the calculations described in 
this paper for Turbospectrum/MARGS. 


[M/H] value. 
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The opacities used in ASSeT are largely independent from those in Kurucz’ model 
atmospheres, especially in the case of continuum (bound-free and free-free) opacity. 

We have, nevertheless, checked that the absolute fluxes predicted in the optical and 
near-infrared are fairly similar for ASSeT and ATLAS9 for a solar-like star, where 
dominates the continuum opacity (as for most APOGEE targets), while in the near-UV the 
differences become larger due to different photoionization cross-sections for metals (mainly 
iron, but also magnesium and other atoms). The overall equation of state is also similar for 
the two codes, since when chemical equilibrium calculations are performed in ASSeT, the 
hnal electron density for solar-like stars is consistent with the original from the ATLAS9 
model at a level of a few percent. In all the ATLAS9 DR12 grids the electron density was 
iterated in ASSeT to be consistent with the atomic and molecular species considered in 
the synthesis. The calculations for a few models with extreme compositions did result in 
iterated electron densities that shifted the tx ~ 100 layer outside of the ATLAS9 structure. 
ASSeT is set to stop when that occurs, so those calculations were successfully completed 
relaxing that limit to ta ~ 10 . 

The ATLAS9 APOGEE GK- and F-class model atmosphere grids mentioned 
above are used in conjunction with the DR12 atomic/molecular line lists (linelist 
201312161124; Shetrone et ah 2015) to compute ATLAS9/ASSeT synthetic spectra with 
Eve microturbulence values (.^ = 0.5, 1, 2, 4, and 8 km s“^), five N abundances ([N/M] 
from —1.0 to -1-1.0 dex in steps of 0.5 dex), and the solar G isotopic ratio (^^G/^^G=89). 
Thus, in the end the GK- and F-class ATLAS9/ASSeT grids contain seven dimensions (7D; 
Teff, log g, [M/H], [G/M], [N/M], and [a/M]) and the same number (1,715,175) of total 
synthetic spectra. The model parameters for the ATLAS9/ASSeT stellar spectral libraries 
are summarized in Table 1. Finally, the synthetic spectra are smoothed to the APOGEE 
resolution, continuum-normalized, re-sampled, and transformed to vacuum wavelengths. 
The ATLAS9/ASSeT spectral libraries are then compressed using Principal Gomponent 
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Analysis method (PCA, Pearson 1901) to prepare it for use by ASPCAP (see Garcia 
Perez et al. 2015 for further details). We note that the line spread function (LSF) for the 
ATLAS9/ASSeT DR12 spectral library is a combo of those derived from several APOGEE 
hbers; actually an average of hve hbers LSFs characterized with a Gauss-Hermite function 
of variable resolution with wavelength (see also Garcia Perez et al. 2015 for further details). 
However, to make a fair comparison with the other spectral libraries treated in this paper 
and to avoid also systematic effects, we are using an ATLAS9/ASSeT spectral library 
version that has been smoothed with a Gaussian kernel to the APOGEE spectral resolution 
{R = 22,500), as the MARGS/Turbospectrum one. 


2.2. MARCS/Turbospectrum spectral library 


The APOGEE MARGS/Turbospectrum spectral library is based on the most 
recent MARGS model atmospheres (Gustafsson et al. 2008). Basically, MARGS model 
atmospheres are one-dimensional models in hydrostatic equilibrium that are calculated 
in LTE and adopting the mixing-length theory for convection (Henyey et al. 1965). The 
MARGS models are spherical for surface gravities log g < 3, while they are plane-parallel 
at higher gravities (appropriate for dwarf stars; see Gustafsson et al. 2008). Line opacities 
are treated with opacity sampling (OS). 


The grid of MARGS model atmospheres for APOGEE was presented in Meszaros et al. 
(2012 )g and we refer the reader to this paper and Gustafsson et al. (2008) for more details. 
MARGS models covering nine Teg values (from 3500 to 5500 K in steps of 250 K)!^, eleven 


®The MARGS model atmospheres can be found in the MARGS Web site; 
http://marcs.astro.uu.se/, 

®We note that the MARGS models with T^g = 3700 K are used instead of those with 
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log g values (from 0 to 5 dex in steps of 0.5 dexfl. seven metallicities [M/H] (from —2.5 to 
+0.5 dex in steps of 0.5 dex), and 25 combinations of C and a-element abundance^lj ([C/M] 
and [a/M] from —1.0 to +1.0 dex in steps of 0.5 dex) are present in the MARCS-APOGEE 
grid used to construct the corresponding MARCS/Turbospectrum spectral library. This 
results in a grid of ~17,325 MARCS models; 1,062 models, however, do not converge and 
are therefore missing in the MARCS-APOGEE grid. These non-converged MARCS model 
atmosphere structures are replaced with the nearest model in chemical space at the same 
Teff and log g (see below). 

Synthetic spectra in the APOGEE spectral range were generated with the 
Turbospectrum package (Alvarez & Plez 1998; Plez 2012), which shares the same input 
data and radiative transfer routines with MARGS (Gustafsson et al. 2008). Turbospectrum 
is a ID LTE spectral synthesis code that accounts for 600 molecules and uses the treatment 
of collisional line broadening described by Anstee Sz O’Mara (1995) and Barklem et al. 
(2000). It allows the computation of flux (or intensity) synthetic spectra for stars of spectral 
type E and cooler. We used version 13.1 of Turbospectrum but modified by us to use van 
der Waals broadening (the Barklem treatment was used when possible, otherwise we used 
van der Waals constants from Kurucz; see Shetrone et al. 2015 for more details) for the 
atomic line^. The synthetic spectra are computed for an array of standard air wavelengths 
with a wavelength step of 0.03 A. By using the MARGS-APOGEE grid (see above) and 


3750 K, because the Tes step in the MARGS grid is 100 K below 4000 K. 

^°Note that plane-parallel and spherical MARGS model atmospheres are computed with 
a microturbulence (.^) of 1 and 2 km s“^, respectively (see, e.g., Meszaros et al. 2012). 
^^The a-elements in MARGS are O, Ne, Mg, Si, S, Ar, Ga, and Ti. 

Version 14.1 of Turbospectrum includes these changes and it is publicly available at 
http://ascl.net/1205.004, 
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the DR12 atomic/molecular line lists (Shetrone et al. 2015) as input in Turbospectrum, we 
constructed MARCS/Turbospectrum synthetic spectra with microturbulent velocities of ^ 
= 0.5, 1, 2, 4, and 8 km s“^ and with varying N content ([N/M] from —1.0 to +1.0 in steps 
of 0.5 dex) and solar ratio. This resulted in a MARCS/Turbospectrum grid with 

seven dimensions (7D; Teg, log g, [M/H], [C/M], [N/M], and [a/M]) containing ~ 43,375 
synthetic spectra; the model parameters for the MARCS/Turbospectrum stellar spectral 
library are also given in Table 1. It is to be noted here that in some cases (~100) the 
spectral synthesis does not converge. However, these mainly correspond to the most extreme 
and unrealistic abundance patterns (e.g., very high or very low a element abundances). As 
in the case of the missing MARCS model atmospheres mentioned above, we also replace 
the missing synthetic spectra (fluxes) with the ones at the same Tes and log g and with the 
nearest chemical composition. Finally, the MARCS/Turbospectrum synthetic spectra are 
processed (i.e., smoothed, continuum-normalized, etc.) and PCA-compressed in the same 
way as the ATLAS9/ASSeT library. 


3. Comparison between ATLAS9/ASSeT — MARCS/Turbospectrum spectral 

libraries 

In this section we compare the Gaussian smoothed GK-ATLAS9/ASSeT and the 
K-MARCS/Turbospectrum spectral libraries, which since these grids overlap in atmospheric 
parameters in the range 3500-5500 K. 

Such comparisons can give indications on the uncertainties due to the use of different 
model atmospheres, different spectral synthesis codes, as well as partly different input data, 
with only the line lists being the same. It also allows to check the possible influence of 
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sphericity effects in the synthetic spectra (Sect. 13.1(1^^ and to explore possible systematic 
differences between both synthetic libraries in the 7D parameter space (Sect. 


A total of 2,552 GK-ATLAS9/ASSeT test synthetic spectra were chosen to uniformly 
span the parameter space, while avoiding the grid boundaries; e.g., by selecting one in ten 
consecutive spectra with different 7D parameters in the nodes of the grid. Our synthetic 
spectra have surface gravities from log g = 0.5 to 4.5 dex, effective temperatures Tes from 
3750 to 5250 K, microturbulent velocities from 1 to 4 km s“^, [C, N, a/M] from —0.5 to 0.5 
dex, and [M/H] from —2.0 to 0.0 dex. 


3.1. Comparisons between spectral syntheses 

For the study of possible differences (e.g.,sphericity effects) between both families 
of synthetic spectra, we extracted the same 2,552 test synthetic spectra from the 
K-MARCS/Turbospectrum spectral library. The differences in the synthetic spectra 
are estimated by deriving the root mean square difference (r.m.s.) between the GK- 
ATLASO/ASSeT and the K-MARGS/Turbospectrum pair of spectral syntheses with the 
same 7D parameters. We performed comparisons for two groups: low-gravity stars (log g < 
2.0) and high-gravity stars (log g > 3.0). In this way, if sphericity effects are affecting the 
H—band synthetic spectra of the low-gravity stars, we would hnd higher r.m.s. values for 
this group. 

Our results are summarized in Figs. [1] and [21 where we compare the r.m.s. values 
from some synthetic spectra in our test sample; i.e., those with log g < 2.0, [G/M] = 0.0, 
[N/M]= 0.0, [a/M] = 0.0 and [M/H] = 0.0 (Fig. [H showing ten synthetic spectra in our 

^^Sphericity effects are expected to be noticeable for low gravity stars, as all ATLAS9 
model atmospheres have plane-parallel geometry. 
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test sample) and those with log g > 3.0, [C/M] = 0.0, [N/M]= 0.0, [a/M] = 0.0 and [M/H] 
= 0.0 (Fig. [21 displaying seven synthetic spectra in onr test sample). The r.m.s. values (of 
the order of ~0.3-0.6 %) are very similar for the two groups. Interestingly, this indicates 
that the overall sphericity effects on the chemical abundances are smaller in the H—band 
than in the optical wavelength region (Heiter & Eriksson 2006 found abundance differences 
as high as 0.35 dex); one reason may be that in the H—band we observe deeper atmospheric 
layers. Indeed, we obtain similar r.m.s. values for other chemical compositions ([C/M], 
[N/M], [a/M] and [M/H] values), validating the adoption of our official DR12 library based 
on ATLAS9 models without any signihcant biases due to the use of plane-parallel model 
atmospheres. The low r.m.s. values obtained are indicative of no large differences between 
the computations based on ASSeT and Turbospectrum spectral synthesis codes. We also 
note that the (systematically) deviating features in Figs. [Hand [2] are hydrogen lines; this 
is because the spectral synthesis codes (ASSeT and Turbospectrum) use different internal 
data for H. However, in the following section the H lines are not used in ASPCAP in htting 
the best £t spectra. 


3.2. Systematic differences 

With the ultimate goal of exploring further the possible systematic differences between 
both grids of APOGEE synthetic spectra, the 2,552 synthetic spectra in our test sample, 
as extracted from the GK-ATLAS9/ASSeT spectral library, have been fitted with the 
K-MARCS/Turbospectrum library using ASPCAP (see Garcia Perez et ah 2015). In 
other words, we find out what MARGS/Turbospectrum seven (7D) stellar parameters are 
recovered by the pipeline when we treat the GK-ATLAS9/ASSeT synthetic spectra as input 
observed spectra. 

The results of this exercise are reported in Table [2] and in Figures 3 and 4. In Table [2] 
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(columns 2 and 3) we show the median in the difference A (= output MARCS parameter 
— input ATLAS9 parameter) and the dispersion (a) of the differences obtained for the 
full sample. To avoid outliers from biasing the statistics, a is computed as the difference 
between the maximum and minimum A after excluding the largest 15.85% of the sample 
and the smallest 15.85%, and divide it by two, which would correspond to the standard 
deviation in a Normal distribution (see Fig. 3). The results of fitting the 2552 test synthetic 
spectra with the GK-ATLAS9/ASSeT library are displayed in Table 2 (columns 4 and 5) 
and Figure 4. In addition, we compare the ATLAS9 versus MARCS 7D parameters for two 
different sub-samples: low-gravity stars (log g < 2.0.) and high-gravity stars (log g > 3.0, 
see Table m columns 6 to 9). 

The MARCS/Turbospectrum spectral library (full sample) provides slightly higher 
effective temperature, microturbulence, [N/M], and [a/M], with median values of 38.1 K, 0.02 
dex, 0.09 dex, and 0.02 dex, respectively. On the other hand, the MARCS/Turbospectrum 
library provides slightly lower metallicity and surface gravity than the ATLAS9/ASSeT 
one; median values of —0.03 dex and —0.13 dex are found for [M/H] and log g, respectively. 
Very similar carbon abundances ([C/M]) are obtained with the MARCS/Turbospectrum 
spectral grid. In order to check if the differences in the 7D parameters mentioned above 
are significant or merely the results of degeneracies, we compare the previous results 
with the ones corresponding to the use of the ATLAS9/ASSeT library as input {observed 
spectra) but also running ASPCAP with the same ATLAS9/ASSeT library. These results 
are reported in Table 2 (cols. 4 and 5), where we a find good consistency for log g and 
Tgff, which have median values of 0.01 ± 0.04 dex and 3.10 ± 17.15 K, respectively. We 
therefore conclude that the systematic differences between MARCS/Turbospectrum and 
ATLAS9/ASSeT for log g and are significant. The nitrogen abundance ([N/M]) may 
be slightly higher (median 0.09) and is the most problematic parameter to recover for the 
MARCS/Turbospectrum library, with the highest a of about 0.20 dex in [N/M] (see Table 
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2). This N problem is not specific to the MARCS/Turbospectrum spectral library, since 
the [N/M] parameter displays also the highest a if we just compare the ATLAS9/ASSeT 
synthetic spectra with themselves (see Table 2 and Fig. 4). We note that ASPCAP is 
limited in accuracy for low-metallicity spectra ([M/H] < —1.0), since the scarcity of lines 
in that regime causes degeneracies among the stellar parameters (see Garcia Perez et al. 
2015). Moreover, the a in N and in the other stellar parameters found in tests by Garcia 
Perez et al. (2015) using ATLASO/ASSeT libraries are similar to the ones found here. 

If we compare the MARGS — ATLAS9 residuals obtained for the 7D parameters by 
surface gravity groups, we hnd that high-gravity stars (log g > 3.0) display higher a values 
(for all parameters) than the low-gravity stars (log g < 2.0, see Table 2; columns 6 to 
9). However, the median parameter values are quite similar with the exception of Teg, 
where the median value for low-gravity stars is signihcantly higher (by about 60 K). The 
differences in the 7D parameters between the two log g groups become evident in Figures 
5 and 6. These hgures show the residuals in surface gravity (A log g) versus the residuals 
in the other stellar parameters (from top to bottom: [G/M], [N/M], [a/M], [M/H], and 
Teff) for low-gravity stars (Fig. 5) and high-gravity stars (Fig. 6). For low-gravity stars 
(Fig. 5), the residuals for most of the data points are close to their median values, with 
the exception of the already mentioned outliers (i.e., the largest and smallest 15.85% of 
the sample) and the N problem. These outliers are dominated by input ATLAS9/ASSeT 
spectra of low-metallicity ([M/H] < —1.0 dex). In spite of the outliers, the results for 
the MARGS/Turbospectrum and ATLAS9/ASSeT spectral libraries in low-gravity stars 
are very similar, and as mentioned above, the influence of possible sphericity effects on 
the derived abundances using ATLAS9 plane-parallel model atmospheres is small in the 
H-band. High-gravity outliers in this analysis correspond mainly to the ATLAS9/ASSeT 
input synthetic spectra with [M/H] < —1.0 dex. Further work by the APOGEE ASPGAP 
team is needed to fully understand why ASPGAP results degrate at high gravities. 
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4. Comparisons of the Solar and Arcturus spectra with Syntheses from 

Spectral libraries 

In order to further investigate how consistent the SDSS-lIl APOGEE spectral 
libraries are, we have htted the Sun and Arcturus observed spectra using the two spectral 
libraries computed with both families of model atmospheres and spectral syntheses codes 
(ATLAS9/ASSeT and MARCS/Turbospectrum). We also compare with the synthetic 
spectra obtained by using ATLAS9 model atmospheres and the MOOgQ spectral synthesis 
code (Sneden 1973). This is relevant because ATLAS9/MOOG synthetic spectra were 
used in the development of the DR12 line lists (see Shetrone et ah 2015). For the Sun, 
in particular, it is of interest to test how good our spectral libraries perform comparing 
with a spectrum at much higher spectral resolution than APOGEE’s, and with well-known 
abundances. In addition, given that APOGEE observes mostly giant stars, we can verify 
how well our synthetic libraries reproduce the molecular lines (i.e., those suitable for GNO 
abundance determinations) using the spectrum of a cooler giant star as Arcturus. 

For the Sun, we use the high-resolution flux spectrum {R = 400,000) by Livingston 
& Wallace (1991). The Sun’s synthetic spectra were computed adopting MARGS and 
ATLAS9 model atmospheres with T^s = 5777 K, log g = 4.4370, solar composition by 
Asplund et ah (2005), and a microturbulent velocity of ^ = 1.1 km s“^. Three different 
synthetic spectra were computed with the following model atmospheres and spectral 
synthesis codes: MARGS/Turbospectrum, ATLAS9/ASSeT, and ATLAS9/MOOG. The 
solar macroturbulent velocity was taken into account in the synthetic spectra by convolving 
them with a Gaussian profile having a FWHM of 1.58 km s“^ (Allende Prieto et ah 2001). 
The synthetic spectra were also convolved with another Gaussian profile to match the 


^^http://www.as.utexas.edu/$\sim$chris/moog.html 
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observed spectrum (FWHM = 1.87 km s“^)0. Finally, the spectra were interpolated to the 
wavelengths (in air) of the observed solar spectrum. All three computed syntheses were 
compared with the solar spectrum. The results from this comparison indicates a fairly good 
agreement between the observed synthetic spectra, as well as a good agreement between 
the synthetic spectra among themselves. The resulting global value for the Sun htting 
were 15.05, 21.43 and 17.70 with the MARCS/Turbospectrum, ATLASO/ASSeT, and 
ATLAS9/MOOG synthetic spectra, respectiveljo. Figure [7] shows the quality of the hts to 
the solar spectrum for the spectral range 16500—16560 A, which includes several atomic 
and molecular lines. The differences (or residuals) between the different types of synthetic 
spectra and the Sun’s observed spectrum are lower than 3% for most data points. The 
r.m.s. value for the MARCS/Turbospectrum synthetic spectrum is slightly lower (~0.1%) 
than the ones for the ATLASO/ASSeT and ATLAS9/MOOG spectra. We note, however, 
that a perfect match between synthetic and observed spectra is not expected because of the 
convective line shifts and asymmetries in the real stars, and of course because our modeling 
of the solar atmosphere is not perfect. 


In the case of the giant star Arcturus, we have used the FTS observed spectrum 
smoothed to the resolution of APOGEE {R = 22,500) and the best ht spectrum from each 
library provided by ASPCAP/FERRE0 (see also next section for further details). The 


^^The macroturbulence of 1.58 km s“^ is from optical spectra (Allende Prieto et ah 2001) 
and we need an extra macroturbulence contribution to match the observed H-band solar 
spectrum. 

^®We estimate the error on the observed spectrum calculating the standard deviation over 
a spectral region free of absorption lines and assuming this error constant along the full 
range in wavelength. 

^^FERRE is available from http: //hebe. as. utexas . edu/f erre, 
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reduced global log value for the fitting was 0.40 for the MARCS/Turbospectrum library 
and 0.41 for the ATLASO/ASSeT one (see Table 3 and Section 5). We focus our comparisons 
of molecular features in those spectral regions selected in the line-by-line abundance analysis 
by Smith et al. (2013). These authors used four windows to extract the abundance 
from 12C160 lines (15578-15586, 15774-15787, 15976-16000, and 16183-16196 A), four 
windows to extract the abundance from lines (15277—15282, 15390—15392, 

15504-15507, and 16189-16193 A), and nine molecular lines of (15260, 15322, 

15397, 15332, 15410, 15447, 15466, 15472, and 15482 A) to extract the abundance of 

The average of the residuals to hts to the observed Arcturus spectrum were derived 
independently in each region/line for the two synthetic libraries. Very small (r.m.s. ~ 

0.1 - 0.3%) differences are found between syntheses from both libraries and the observed 
spectrum (and between the library synthetic spectra themselves). Figures IS and [9] display 
the hts in the and spectral regions and, in the spectral lines mentioned 

above. In the windows, the MARCS/Turbospectrum library spectrum hts slightly 

better the regions at 15578—15586 and 15976—16000 A, while the ATLASO/ASSeT 
library spectrum hts better the regions at 15774—15787 and 16183 — 16196 A. For the 
windows, the MARCS/Turbospectrum library spectrum hts better the regions 
at 15277—15282 and 15390—15392 A, while the ATLASO/ASSeT library spectrum hts 
better the regions at 15504—15507 and 16189—16193 A. Finally, for the lines, the 

MARCS/Turbospectrum library spectrum hts slightly better the lines at 15260, 15397, 
15466, and 15482 A, while the ATLASO/ASSeT library spectrum hts better the lines at 
15322, 15332, 15447, and 15472 A. Both libraries provide just the same residual average 
for the CN line at 15410 A. We conclude that both synthetic spectral libraries provide 
an excellent ht to the molecular features in the spectrum of Arcturus; the r.m.s. values 
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between both synthetic libraries are no signihcant and of the order of only ~ 0.1 



5. Application of the APOGEE spectral libraries to selected giant stars 

Given that most of the APOGEE sample are red giant stars, we are interested here 
in exploring the abnndance differences, as obtained by the APOGEE spectral libraries, in 
a small sample of well known giant stars observed at very high resolntion. Smith et ah 
(2013) derived chemical abnndances from a line-by-line analysis of 15 elements in several 
well-known bright held giants and explored what elements can be analyzed from APOGEE 
spectra. The sample analyzed here consists in the fonr Smith et ah (2013) stars with 
> 3500 K; this inclndes two M-giants {(3 And and S Oph) and two K-giants (a Boo and 
fi Leo). For their stndy. Smith et ah nsed high-r esolution spectr a in the H-band acquired 


with the Fourier transform spectrometer (FTS; 


Hall et ah 


19791) installed at the Goude 


focus of the Kitt Peak National Observatory 4m Mayall rehector. The spectral resolution of 
these FTS spectra varies from 45,000 to 100,000. Th e origi nal spectra cover a wavelength 


range larger than that of APOGEE, but 
sp ectral range 


Smith et ah 


(120131) restricted their analysis to the 


by 


Hinkle et al 


r om 15 00 to 1700 nm. The infrared atlas spectrum of a Boo (Arcturus) 
fjl995l) . obtained with the same instrument at a resolution of 100,000, is 


added to our sample of stars with FTS spectra. The data were smoothed to the APOGEE 
resolution (i.e., R = 22,500) by convolving with a Gaussian kernel. We also converted the 
resulting convolved spectra to the APOGEE apStar FITS format described by Holtzman et 
al. (2015). Our FTS sample includes the H—band spectra of a Boo and fi Leo obtained 
with the APOGEE spectrograph, but using the New Mexico State University 1.0-meter 


^®Note that the f-values for atomic lines were tuned to match the Sun and Arcturus but 


not those of molecular transitions (see Shetrone et al. 2015). 
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Telescope (Holtzman et al. 2015; NMSU Im, hereafter). The NMSU Im spectra are 
reduced with the APOGEE data reduction pipeline (Nidever et al. 2015). All spectra 
mentioned above were processed with a qnick version of ASPCAP, QASPCAP, which is 
short version of ASPCAP, and prepares the spectra for the antomated htting with EERRE 
(see Garcia Perez et al. 2015 for a detailed description of ASPCAP). Finally, we derived 
the atmospheric parameters and chemical abnndances of 15 elements (see below) with 
FERRE, interpolating in the ATLASO/ASSeT and MARCS/Turbospectrum synthetic grids 
described in this paper. 

The atmospheric parameters and chemical abnndances obtained are listed in Tables 0] 
and[5l respectively. For all stars, we hnd a very good agreement between the valnes obtained 
with both synthetic grids. The log of th® fits for each spectral library are qnite similar 
(Table [3]). The differences of the derived Teg, log g, microtnrbulent velocity (in log scale), 
and [M/H] for each star (and synthetic grid) are plotted in Figure [TOl The atmospheric 
parameters derived here for the Arcturus giant are somewhat different from those adopted 
in Smith et al. (2013): ATeff < 90 K, Alog g < 0.4 dex, A^ < 0.6 km s“^, A[M/H] < 0.2 
dex. In Smith et al. (2013) the effective temperatnres were based on the (J—K) color and 
derived from an average of two calibrations: Gonzalez Hernandez & Bonifacio (2009) and 
Bessell et al. (1998). Their snrface gravities were obtained from evolutionary tracks and 
microtnrbulent velocities from the Fe I lines. The stellar parameters derived here are pnrely 
spectroscopic. In addition, it is important to note that Smith et ah (2013) carried out a 
line-by-line chemical abundance analysis using the MOOG synthesis code and a so-called 
intermediate version of the APOGEE line list (line list INT] Shetrone et al. 2015) that is 
previous to the DR12 APOGEE line list used here. 

Regarding abnndances, we hnd in general very good agreement (<0.1 dex) between 
the chemical abnndances obtained by the ATLASO/ASSeT and MARGS/Tnrbospectrnm 
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stellar spectral libraries in the FTS stars (see Fig. [n]). For Arcturus {a Boo), 
MARCS/Turbospectrum performs slightly better than ATLAS9/ASSeT, with the exception 
of N (see below). Despite the higher MARCS/Turbospectrum Tgff values, the set of derived 
MARCS/Turbospectrum abundances in Arcturus (and also in /i Leo) are slightly lower 
(<0.1 dex) than those from ATLAS9/ASSeT because of the generally lower metallicity and 
gravity obtained (which compensates the expected abundance increase due to a higher 
Tgfj; Section 3) with the MARCS/Turbospectrum library (see Tables 0] and |5]) . However, 

N seems to be more affected and the MARCS/Turbospectrum derived N abundances 
can be lower by 0.15 dex. On the other hand, ATLAS9/ASSeT hts slightly better than 
MARCS/Turbospectrum in the two FTS stars (3 And and 6 Oph; especially concerning 
nitrogen, where MARCS/Turbospectrum gives N abundances lower by ~0.2—0.3 dex (see 
Fig. [TT]) . This, however, is likely due to the use of the cooler 3700 K MARCS model in 
place of the 3750 K one, which is lacking in the MARCS grid, to interpolate their exact 
matching value of effective temperatures of about 3825-3850 K. 


In short, based on the results from the comparisons performed in this study, there 
is good indication that the MARCS/Turbospectrum library, although with several more 
models missing from the grid, gives results comparable to those from the ATLAS9/ASSeT 
library 1^. 


6. Conclusions and future work 

We present the stellar spectral libraries for the hnal data release of the SDSS-III 
APOGEE survey, which are used for the automated chemical analysis of survey data. The 

^^Garcia Perez et ah (2015) used the same ATLAS9/ASSeT spectral library to analyze 
the spectra of FTS stars but with a different order in the stellar parameters. 
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spectral libraries employed in the data release 12 (DR12; Alam et al. 2015) are constructed 
for a wide range in effective temperature (Tgff ranging from 3500 to 8000 K) and are based 
on ATLAS9 model atmospheres and the ASSeT spectral synthesis code. We also present 
here a second family of SDSS-III APOGEE stellar spectral libraries based on MARCS 
model atmospheres and the Turbospectrum spectral synthesis code. The ATLASO/ASSeT 
(Teff = 3500-8000 K) and MARCS/Turbospectrum (Teg = 3500-5500 K) synthetic grids 
have seven dimensions (7D), covering a wide metallicity ([M/H]), surface gravity (log g), 
microturbulence (.^), carbon ([C/M]), nitrogen ([N/M]), and a-element ([a/M]) ranges of 
variation. 

We have compared both ATLASO/ASSeT and MARCS/Turbospectrum spectral 
libraries to a test sample of 2552 synthetic spectra with the same 7D stellar parameters. 
The differences found between both families of synthetic spectra are very small (r.m.s. 
values of the order of only ~0.3—0.6 %). Interestingly, we hnd that the sphericity effects 
in the H—band seem to be smaller than those previously found in the optical range and 
the ASSeT and Turbospectrum spectral synthesis codes provide very similar synthetic 
spectra. By fitting the GK-ATLASO/ASSeT library with the K-MARCS/Turbospectrum 
library, we have found small systematic differences in the seven main stellar parameters 
(7D; Teff, [M/H], log g, [C/M], [N/M], and [ct/M]) automatically provided by the 
SDSS-III APOGEE survey for low-gravity stars (log g < 2.0). The outliers correspond to 
low-metallicity ([M/H] < —1.0) synthetic spectra. However, the results for high-gravity 
(log g > 3.0) stars are worse than the former ones, and the average scatter for the entire 
parameter space is higher than in low-gravity stars. These outliers are also dominated by 
synthetic spectra with [M/H] < —1.0 dex. Further work by the APOGEE ASPCAP team is 
needed to completely understand the presence of these outliers in the 7D parameter space. 


Both the DR12 SDSS-III APOGEE synthetic spectral library as well as the additional 
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spectral library based on the MARCS model atmospheres provide almost identical model 
fits to the observed spectra of the Snn, Arctnrns, and the stars with FTS spectra. For 
example, they give an excellent fit to the Snn’s spectrnm as well as to the molecnlar 
featnres (CO, OH, and CN) in the spectrnm of Arctnrns; the differences (or residuals) 
between these synthetic libraries are of the order of only ~ 0.1 - 0.3% (r.m.s.). We conclude 
that both SDSS-III APOGEE synthetic spectral libraries provide very similar results 
(i.e., atmospheric parameters and chemical abundances), which supports the use of the 
ATLASO/ASSeT synthetic grids (which otherwise cover a parameter space much wider 
than the actual MARCS/Turbospectrum grid) in DR12. The SDSS-III APOGEE synthetic 
spectral libraries presented here are publicly available online and they can be used also for 
chemical analysis in the H—band making use of other available high-resolution spectroscopic 
instruments working in the H—band. 

The APOGEE stellar spectral libraries presented here will be improved for the 
SDSS-IV/APOGEE-2 survey, and periodically updated in the future. We plan to extend 
the MARGS/Turbospectrum stellar spectral library to cooler temperatures (2500 < < 

3500 K). The effect of other molecules such as H 2 O and FeH may be important at these 
extremely cool effective temperatures and we will need to update the present APOGEE 
linelist by including these molecules. Finally, we plan to evaluate the effects of the missing 
opacities for polyatomic molecules (like HGN, G 2 H 2 ) on the structures of the cool (Teg < 
4000 K) and G-rich MARGS models atmospheres and we plan to improve the latter G-rich 
models with new opacities for such polyatomic molecules. 

The authors acknowledge the Texas Advanced Gomputing Genter (TAGG) at The Uni¬ 
versity of Texas at Austin for providing HPG resources that have contributed to the research 
results reported within this paper (URL: http://www.tacc.utexas.edu). This paper made use 
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Table 1. Parameters of the ATLAS9/ASSeT and MARCS/Turbospectrum spectral libraries 


Class 

Teff range (step) 

log g range (step) 

[M/H] range (step) 

[C/M] range (step) 

[N/M] range (step) 

[a/M] range (step) 

log ^ range (step)^ 


(K) 

(dex) 

(dex) 

(dex) 

(dex) 

(dex) 

(km s“^) 

ATLAS9/ASSeT 







1 

GK 

3500-6000 (250) 

0-5 (0.5) 

-2.5-0.5 (0.5) 

-1.0-1.0 (0.25) 

-1.0-1.0 (0.5) 

-1.0-1.0 (0.25) 

-0.301-0.903 (0.30l^ 

F 

5500-8000 (250) 

0-5 (0.5) 

-2.5-0.5 (0.5) 

-1.0-1.0 (0.25) 

-1.0-1.0 (0.5) 

-1.0-1.0 (0.25) 

-0.301-0.903 (0.301) 

MARCS/Turbospectrum 

K 

3500-5500 (250) 

0-5 (0.5) 

-2.5-0.5 (0.5) 

-1.0-1.0 (0.5) 

-1.0-1.0 (0.5) 

-1.0-1.0 (0.5) 

-0.301-0.903 (0.301) 


The microturbulence (log step is given in logio units (uniform step of 0.301 dex). 












ATLAS9/ASSET vs. MARCS/Turbospectrum synthetic spectr 
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Fig. 1.— Comparison between synthetic spectra extracted from the GK-ATLAS9/ASSeT and K- 

MARCS/Turbospectrum libraries for low-gravity stars (log g < 2). The wavelength ranges covered by the three 
APOGEE detectors are showed. Only ATLAS9/ASSeT spectra are displayed (in black). The residuals, computed as 
MARGS — ATLAS9 fluxes, have been multiplied by a factor Eve of to make them visible in the hgure (red line). The 
7D parameters of each spectra are indicated above each spectrum (where ^ = Vdop), together with the root mean square 
(r.m.s.) value computed from each pair of synthetic spectra with the same 7D parameters. 






































Table 2. Systematic differences between ATLAS9/ASSeT (input) and 
MARCS/Turbospectrum (output) spectral syntheses in the 7D parameter space 


Parameter 

Median(MARCS-ATLAS9) 

Full sample 

(T 

Median(ATLAS9-ATLAS9)* 

Full sample 

(7 

Median(MARCS-ATLAS9) 

log g < 2.0 subsample 

(7 

Median(MARCS-ATLAS9) 

log g > 3.0 subsample 

(7 

log ? 

0.02 

0.08 

0.00 

0.03 

0.01 

0.02 

0.08 

0.19 

[C/M] 

0.00 

0.08 

0.00 

0.04 

0.00 

0.07 

-0.01 

0.10 

[N/M] 

0.09 

0.20 

0.01 

0.15 

0.10 

0.17 

0.07 

0.27 

[ck/M] 

0.02 

0.02 

0.00 

0.01 

0.03 

0.02 

0.01 

0.03 

[M/H] 

-0.03 

0.04 

0.00 

0.01 

0.01 

0.04 

-0.04 

0.04 

log g 

-0.13 

0.15 

0.01 

0.04 

-0.11 

0.13 

-0.11 

0.14 

Tett 

38.10 

50.15 

3.10 

17.15 

64.30 

41.70 

20.00 

59.20 


*These are the results using synthetic spectra from the ATLASQ/ASSeT library as input and output for ASPCAP. 


Table 3. ATLAS9/ASSeT vs. MARCS/Turbospectrum log values in FTS stars 


Star 

log ATLAS9/ASSeT 

log MARCS/Turbospectrum 

a Boo (FTS) 

0.4075 

0.3972 

a Boo (atlas) 

0.3524 

0.3443 

a Boo (NMSU Im) 

2.0242 

2.0236 

IM Leo (FTS) 

1.1852 

1.1885 

fM Leo (NMSU Im) 

1.8115 

1.8260 

/3 And (FTS)'* 

1.3842 

1.3926 

5 Oph (FTS)** 

1.1643 

1.1730 


^These stars have atmospheric parameters corresponding to a hole in the 
MARCS/Turbospectrum grid. 
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Table 4. ATLAS9/ASSeT vs. MARCS/Turbospectrum stellar parameters in FTS stars 



ATLAS9/ASSeT 


MARCS/Turbospectrum 


Star 

Teff 

log g 

[M/H] 


TeS 

log g 

[M/H] 

? 

a Boo (FTS) 

4187 

2.04 

-0.40 

2.03 

4192 

1.85 

-0.47 

2.07 

a Boo (atlas) 

4188 

2.07 

-0.43 

1.90 

4192 

1.95 

-0.48 

1.92 

a Boo (NMSU Im) 

4208 

2.07 

-0.50 

1.15 

4223 

1.92 

-0.56 

1.25 

fi Leo (FTS) 

4493 

2.80 

0.44 

1.93 

4520 

2.76 

0.40 

1.99 

fi Leo (NMSU Im) 

4560 

2.98 

0.36 

1.00 

4551 

2.87 

0.30 

1.04 

13 And (FTS)=- 

3823 

1.16 

-0.20 

2.36 

3791 

1.19 

-0.24 

2.42 

5 Oph (FTS)^ 

3832 

1.45 

0.03 

2.21 

3809 

1.48 

-0.04 

2.31 


^These stars have atmospheric parameters corresponding to a hole in the MARCS/Turbospectrum grid. 









Table 5. ATLAS9/ASSeT vs. MARCS/Turbospectrum element abundances in FTS stars 


Star 

[Fe/H] 

[C/H] 

[N/H] 

[O/H] 

[Mg/H] 

[Al/H] 

[Si/H] 

[K/H] 

[Ca/H] 

[Ti/H] 

[V/H] 

[Mn/H] 

[Ni/H] 

ATLAS9/ASSeT 

a Boo (FTS) 

-0.42 

-0.29 

-0.30 

-0.16 

-0.23 

-0.24 

0.01 

-0.44 

-0.37 

-0.43 

-0.53 

-0.48 

-0.38 

a Boo (atlas) 

-0.45 

-0.31 

-0.35 

-0.20 

-0.25 

-0.37 

-0.10 

-0.45 

-0.42 

-0.43 

-0.58 

-0.49 

-0.39 

a Boo (NMSU Im) 

-0.52 

-0.46 

-0.58 

-0.39 

-0.39 

-0.27 

-0.30 

-0.51 

-0.51 

-0.26 

-0.71 

-0.55 

-0.32 

IM Leo (FTS) 

0.44 

0.43 

0.91 

0.49 

0.36 

0.50 

0.53 

0.50 

0.23 

0.62 

0.30 

0.50 

0.50 

IM Leo (NMSU Im) 

0.31 

0.37 

0.73 

0.37 

0.22 

0.41 

0.33 

0.19 

0.21 

0.64 

0.32 

0.50 

0.^7 

/3 And (FTS) 

-0.21 

-0.35 

0.20 

-0.11 

-0.07 

-0.12 

-0.06 

-0.30 

-0.26 

-0.07 

-0.23 

-0.11 

-0^ 

<5 Oph (FTS) 

0.00 

-0.08 

0.30 

0.08 

0.06 

0.19 

0.24 

0.27 

-0.08 

0.17 

-0.04 

0.21 

o.ds 

MARCS/Turbospectrum 

a Boo (FTS) 

-0.49 

-0.37 

-0.38 

-0.17 

-0.26 

-0.24 

-0.06 

-0.45 

-0.37 

-0.45 

-0.56 

-0.53 

-0.43 

a Boo (atlas) 

-0.51 

-0.36 

-0.50 

-0.22 

-0.28 

-0.37 

-0.15 

-0.47 

-0.42 

-0.43 

-0.59 

-0.53 

-0.42 

«Boo (NMSU Im) 

-0.59 

-0.50 

-0.64 

-0.39 

-0.40 

-0.26 

-0.36 

-0.52 

-0.51 

-0.35 

-0.72 

-0.58 

-0.37 

IM Leo (FTS) 

0.38 

0.40 

0.91 

0.47 

0.32 

0.50 

0.41 

0.50 

0.23 

0.58 

0.30 

0.50 

0.50 

IM Leo (NMSU Im) 

0.24 

0.31 

0.73 

0.32 

0.20 

0.41 

0.27 

0.15 

0.18 

0.57 

0.26 

0.50 

0.41 

/3 And (FTS)'‘ 

-0.27 

-0.44 

-0.08 

-0.34 

-0.03 

-0.29 

-0.02 

-0.36 

-0.38 

-0.28 

-0.36 

-0.15 

-0.23 

(5 Oph (FTS)=' 

-0.07 

-0.16 

0.07 

-0.11 

-0.11 

-0.02 

0.27 

0.19 

-0.19 

-0.08 

-0.15 

0.17 

-0.02 


These stars have atmospheric parameters corresponding to a hole in the MARCS/Turbospectrum grid. 
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Fig. 3.— Systematic differences between the ATLAS9/ASSeT and the 

MARCS/Turbospectrum spectral libraries in the 7D parameter space. The input 
stellar parameters are those from the ATLAS9/ASSeT spectral library, while the output 
parameters are those derived/recovered with the MARCS/Turbospectrum library by using 
ASPCAP. The distribution of the differences (i.e., output MARCS parameter — input 
ATLAS9 parameter) is shown for the 7D grid (in black). A Gaussian distribution £t, 
excluding the ~15.85 % of the largest data differences (i.e., outliers), is also displayed (grey 
curve). 
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Fig. 4.— Same as Figure 3 but comparing both the input and the output from the 
ATLAS9/ASSeT spectral library in the 7D parameter space. 
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A logg (MARCS - ATLAS9) vs. A porometer (MARCS - ATLAS9) in low-grovity stors 
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Fig. 5.— Parameter differences obtained between MARCS/Turbospectrum and 

ATLAS9/ASSeT spectral libraries for low-gravity (log g < 2.0) stars. From top to bot¬ 
tom: differences (MARCS — ATLAS9) in carbon ([C/M]), nitrogen ([N/M]), a—elements 
([a/M]), metallicity ([M/H]), microturbulent velocity (^), and effective temperature (Teff) 
versus MARCS — ATLAS9 differences in surface gravity (A log g). The red line indicates 
the median values of the parameter differences. 
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A logg (MARCS - ATLAS9) vs, A parameter (MARCS - ATLAS9) in high-gravlty stars 
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Fig. 6.— Same as Figure 5 but for high-gravity stars; i.e., log g > 3.0. 
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1.650x10^ 1.652x10^ 1.654x10^ 1.656x10^ 

Wavelength (Angstroms) 

Fig. 7.— High-resolution H—band observed spectrum of the Sun (black dots) in the 
16500—16560 A region versus the best hts obtained using MARCS/Turbospectrum (red line), 
ATLASO/ASSeT (blue line), and ATLAS9/MOOG (green line) synthetic spectra. All spectra 
are expressed in air wavelengths. The residuals, computed as flux(synthetic—observed)-|-0.45, 
are plotted at the bottom with the same colors. The spectral features identihed by Hinkle 
et ah (1995) are indicated at the top. 
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Fig. 8.— Arcturus FTS observed spectrum (black dots) smoothed to the APOGEE resolution 
{R = 22,500), showing four molecular windows/regions. The best hts obtained with 

the MARCS/Turbospectrum (red line) and ATLASO/ASSeT (blue line) synthetic spectra 
are also shown. All spectra are expressed in air wavelengths. The residuals, computed as 
flux(synthetic—observed)+0.60, are plotted at the bottom with the same color code. The 
spectral features identihed by Hinkle et ah (1995) are indicated at the top. 
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Fig. 9.— Arcturus FTS observed spectrum (black dots) smoothed to the APOGEE resolution 
{R = 22,500) and the synthetic spectra in four molecular windows/regions (upper 

panel) and in nine lines (lower panel). Symbols and colors as in Figure 8. The 

spectral features identihed by Hinkle et ah (1995) are indicated at the top. 

















































































44 


a Boo (FTS) a Boo (atlas) 


0.350 

^ 0.175 


0.350 

^ 0.175 


3 

rlj 0.000 


3 

rlj 0.000 


1$ -0.175 


< -0.175 


-0.350 

[M/H] = -0.47; Telf=4275; log g = 1-70; Vdop=1-85 

-0.350 

[M/H] = -0.47; Telf=4275; log g=1-70; Vdop=1-85 


[Mj?H] Teff/100 logg Vdop Teff/100 logg Vdop 


0.350 

0-175 

0.000 

-0.175 

-0-350 

Teff/100 logg Vdop 

IX Leo (NMSU 1 m) 




[M/h] Teff/100 logg Vdop 

0.350 

0-175 

0.000 

-0.175 

-0.350 

[M/H] Teff/100 logg Vdop 



jx Leo (FTS) 



And (FTS) 



Fig. 10.— Comparison between the atmospheric parameters derived by FERRE using 
ATLASO/ASSeT and MARCS/Turbospectrum for the calibration sample observed with the 
FTS. In order to have all the values in the same scale, the value has been divided by 100. 
Moreover, the microturbulent velocity (.^ = 'v^op) is in logarithmic scale, while the dashed 
line is the zero point. Note that stars {3 And and 5 Oph have atmospheric parameters that 
correspond to a hole in the MARCS/Turbospectrum grid (see text). 
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Fig. 11.— Same as Fig. [10] but for the abundances of individual elements. Note that 
stars (3 And and 5 Oph have atmospheric parameters that correspond to a hole in the 


MARCS/Turbospectrum grid (see text). 



















